Bloch electron in a magnetic field and the Ising model 
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The spectral determinant det{II ~el) of the Azbel-Hofstadter Hamiltonian H is related to Onsager's 
partition function of the 2D Ising model for any value of magnetic flux "1> = 2nP/Q through an 
elementary cell, where P and Q are coprime integers. The band edges of H correspond to the 
critical temperature of the Ising model; the spectral determinant at these (and other points defined 
in a certain similar way) is independent of P. A connection of the mean of Lyapunov exponents to 
the asymptotic (large Q) bandwidth is indicated. 



Although the problems of an electron in a constant 
magnetic field and that of an electron on a periodic lat- 
tice were solved in the early days of quantum mechan- 
ics, the case where a magnetic field and a lattice are 
present simultaneously still defies adequate understand- 
ing. The simplest model for this case is commonly re- 
ferred to as the Hofstadter or Azbel-Hofstadter model 
10,^. The corresponding Hamiltonian describes an elec- 
tron on a two-dimensional (2D) N x N square lattice 
with nearest-neighbour hopping subject to a perpendic- 
ular uniform magnetic field: 



nx,ny = Q,l,...,N-l, 



(1) 



where $ = 2t:P/Q is the flux through an elementary 
cell (measured in units of the elementary flux), P and Q 
are coprime (that is they do not have a common divisor 
other than 1) positive integers, A = ti/i2 > is a ra- 
tio of hopping amplitudes in the x and y directions. We 
impose the periodic boundary conditions {ipN,ny — ^o.uyi 
Tp-i^ny — V'Af-i,Tij,j similarly for Uy ). Usually, the Hamil- 
tonian (^) is defined on an infinite lattice at the outset; 
in the present paper, however, we shall deal first with a 
finite N and take the limit — > 00 later on. 

This Hamiltonian and the related one-dimensional op- 
erator were studied in many works (see [^-^ for reviews) . 
Nevertheless, the quantitative description of the spec- 
trum remains to a large extent an open problem. Hamil- 
tonian (1) is related to that of a quantum particle on a 
line in the presence of two periodic potentials. It is the 
effect of incommensurability of the periods (when $/27r 
is irrational or, in other words, P,Q ^ 00) leading to the 
fractal structure of the spectrum that makes the problem 
both difficult and interesting. Its solution would be an 
essential contribution to the theory of fractals, would be 
important for studies of localization-delocalization phe- 
nomena, the quantum Hall effect, as well as purely for 
number theory and functional analysis. 



It was noticed by Wiegmann and Zabrodin Q that 
the model is related to integrable systems and the algebra 
Uq{sl2). This was generalized in In the present paper, 
we find a different type of relation to integrable systems. 
We shall see that the spectral determinant Aet{H — el) 
for any $ = 2i:P/Q is simply mapped onto Onsager's 
partition function of the 2D Ising model (formulas (|^)- 
(pO|)). This relation is, in a sense, complementary to the 
one established in : we shall explain this later in the 
text. 

In the present paper, we remove part of the mystery 
about the Thouless conjecture on the total bandwidth. 
It is well known that the spectrum of H for an infinite 
lattice (on hiZ)) consists of Q intervals (bands). When 
Q ^ 00, the total width W of the bands for A ^ 1 is 
known to approach 4|1 — A| exponentially fast in Q; 
for A = 1, it is of order 1/Q Thouless formulated 
a conjecture that for A = 1 the total width W ^ 

32G/7rQ as Q cx), where G = 1-1/3^ + 1/5^ -1/7^ .. . 
is Catalan's constant. (Notation A ^ B means, hence- 
forth, that A/B tends to 1 in the limit.) This result was 
derived only for P = 1, P = 2 so far jrT|-|l^ but is sup- 
posed to hold in the general case (even when P/ Q tends 
to any nonzero limit) for which it is supported by exten- 
sive numerical data; thus, the bandwidth is supposed to 
have an interesting universality property with respect to 
We shall generalize this conjecture and recast it into 
a simpler type of statement (expression (p^)) by relat- 
ing the quantity 7(e) — fimAr^oo(l/^^) I det(_ff — e/)| 
(which is interpreted as a mean of Lyapunov exponents) 
to the bandwidth. In particular, for any band edge Se^, 
any coprime P, Q, and A = 1, we get an exact result: 
7(^ei) = '^G/ttQ. This differs from the supposed asymp- 
totic formula for the total bandwidth only by the factor 8. 

First, let us use the translational symmetry of H in the 



X and y directions. Substitute ipn^,ny 
2Tik/N, fc = 0, 1, . . . , A^-1 into (1) to get H 
the eigenvalue equation for being 



N-l rr 

fe=0 ^fe' 



1 



f 2'kP 2'Kk\ 

n-l + ^J^n+l + 2AC0S ( -g-"- + ) ^" ^ 



0,1, 



,iV-l. 



When n is allowed to range from — oo to oo, the corre- 
sponding Hk on 12{Z) is called the almost Mathieu oper- 
ator. 

Let us assume that N is divisible by Q and substitute 



2»-'e,, J = 0,1,...,Q- 1, io = ^ 
0, 1, . . . , iV/Q - 1 into (I). We get iJ^ = (S)'^^^^ 



m - 



{HkmOj = + + 2Acos 
j = 0,l,...,Q-l; eQ = e^"6 



27rfc 



(3) 



The Chambers formula ||l5[ gives the dependence of 
the spectral determinant (characteristic polynomial) of 
Hkm on k and m, namely: 



det{Hkm — sIqxq) — 
(^(T(e) - 2A'3 cos 



27rfc 



27rTO \ 



(4) 



where cr(e) is a polynomial of degree Q which depends 
on A and but not on k and m. It is easy to verify (Q) 
by considering the matrix given by (j^) and the one ob- 
tained after the substitution = Yl2^o^ ^^P 2ni{jlP/Q + 
jm/N + lk/N)£^i. It is the expression (^) which im- 
plies that the spectrum of H in the limit — s- oo 
consists exactly of Q bands: the image of the interval 
[-2(1 + A'3),2(l + XQ)] under the inverse of the trans- 
form (T = cr(e) (see Fig. m. 









1 cT(e) 





















FIG. 1. The polynomial a{e) is sketched for Q = 3. The 
bands of the spectrum are shown by thick lines. 



Thus, 

N~lN/Q-l 

dct(i/ - £/) = n n det(i/fc„ - £/qxq) = 

fc=0 m=0 

(-1)"^ (-(^) - 2A^ cos IS. - 2 cos ^ " 



N/Q 



n/q) 



(5) 



In the limit of infinite lattice {N oo), let us replace 
2nk/N in (H) by a continuous parameter 0, denote by 
iJe, and consider the mean 7(e) of Lyapunov exponents 
over all real 9. By virtue of the 3-term recursion (||), 
any /i„ is obtained from initial conditions /io, /ii. The 
Lyapunov exponent 79(2) corresponding to Hg describes 
the exponential rate of growth (or decay) of a solution to 
the equation {Hoii)n — £Mni n = 0, 1, . . .. More precisely. 



70(e) 



lim n 



max 



^l/2 



(6) 



Note that the difference between Hg defined in the space 
of square-summable sequences with indices n = 0, 1, . . . 
and the almost Mathieu operator corresponding to n = 
. . . , — 1, 0, 1, . . . is that the spectrum of the latter is dou- 
bly degenerate. However, the normalized to unity density 
of states pe{x) is the same in both cases. According to 
the Thouless formula |16l , 79 (e) — ^ hi\e — x\pg{x)dx. 
Using this formula and (^, we get for the mean of 76/ (e) 
over 9 (to ensure that the determinant is nonzero we as- 
sume that e has an imaginary part; after taking the limit 
iV ^ 00 we can let Sje ^ 0): 



7(e) = lim i^ln|det(iI-£/)| = 

^5— / / In |o-(£) — 2A'^ cosa; — 2cosy|dxdy. 



(7) 



On the other hand, Onsager's partition function Z for 
the 2D Ising model on a square N x N lattice satisfies 
IitI: 



lim — - In Z ln(2 sinh a') — 



1 

2^ 



"'0 



In 



smh a 
2-— — 7 cos a; — 2 cosy 
smh a 



^ cosh a cosh a' 
sinh a' 

dxdy. 



(8) 



Here a = 2J/T, a' = 2J'/T, where T is the tem- 
perature; J, J', interaction constants in x and y di- 
rection. We now set sinh a/ sinh a' = A'^ (hence a' = 
arcsinh (A^'^ sinh a)); 



a{e) = 2 



cosh a cosh a' 
sinh a' 



^ 2A^ cothaVl + A-2Q sinh^ a; 



(9) 



2 



to obtain the equivalence (as N 
coprime P and Q: 



oo) for any e and 



\ det t2{H ~el)\ - Z^/Q, 



(10) 



where ^ 2 sinh a'. 

A simple analysis of (||) shows that for any real T, 
|c(£)| > 2(1 + A'^). The minimum of |cr(e)| as a func- 
tion of T, \(j{e)\ = 2(1 + A'^), corresponds to the critical 
temperature Tc of the Ising model (at Tc the argument 
of the logarithm in (|^) and (H) vanishes at an integra- 
tion limit). On the other hand, the energies Eg. at which 
|fT(£ei)| = 2(1 -t- A'^) are the band edges of limAr_+oo H. 

It is well known |l^ that Z at Tc equals the square 
of the partition function Zjj of dimers on the lattice if 
ri = V 2 sinh a and r2 — V 2 sinh a' are interpreted as 
activities of dimers in x and y directions. Let us identify 



1 — 



.2 _ ,Q 



(recall that ti and t2 are hopping 
amplitudes in x and y directions, and A = ^1/^2)- One 
result of Lieb and Loss |l9j says that det H ^ Z\, as 
— > 00 for the flux $ = tt, that is for P = 1, Q = 2. 
Here H ~ t2H — a more symmetric form of the Azbel- 
Hofstadter Hamiltonian. Since it is known that e = is a 
band edge for Q even, we have from (^^ det H ~ Z^^ , 
which is a generalization of the mentioned result to any 
coprime P odd, Q even. 

In the works the Bethe ansatz equations are for- 
mulated for the zeros of the characteristic polynomial 
of Hkrm in other words, for the roots of the equation 
a{e) = const. Thus are to do with the "internal" 
structure of cr(e), whereas in our formalism <7{e) enters as 
an "indivisible" object and our results come from what 
appears (at least at first sight) "external" to a{e) struc- 
ture of the problem: the translational symmetry and the 
Chambers formula. 

Henceforth, we only consider the limit of infinite A^. 
Let us simplify expression (^ (note an important fact 
that 7 depends on e and P only through a{e) !). Because 
of the obvious equality 7(0-, A) = In A + ^{a/X'^ , V'^); it 
is sufficient to consider only the case < A < 1. In this 
case, taking the integral over y in (|^) we get: 

7(cr) = 

'' Trb /o^ 8'rccosh(|a-|/2 -1- A*^ cosa;)da;, 
|a|>2(l + AQ); 
j_ ^™{(2-|.|)/2A«} ^^^^^^^ + \Q cosx)dx, 



2(1 + AQ) > \a\ > 2(1- A'?); 
[0, |a| <2(1-AQ). 



(11) 



The fact that 7(e) is zero on the image of the interval 
a e [-2(1 - A'3),2(l - A'5)] is in agreement with the 
general argument: this image is an intersection of the 
intervals of the spectra of the operators Hg. The gen- 
eralized eigenfunction of any Hg is just a Bloch wave 



on a periodic lattice with Q atoms in an elementary cell. 
Therefore, the Lyapunov exponents "fg (e) are zero on this 
set. By definition, 7(e) is a mean of these exponents. 

Let us remark that rather than considering det{H — 
£/), it is possible to derive equations (11) making use of 
some known properties (see, e.g., [Q) of the monodromy 
operator. 

Now we can note that in the limit Q —^ 00, a fixed, 
7(ct) tends to the following Lyapunov exponent on the 
spectrum of H: 7(e) ^ if A < 1, 7(e) — InA if A > 1, 
in accordance with a statement of Aubry and Andre [^ . 
(Naturally, we obtain the same asymptotic result for any 
individual ^g{e).) 

For Tc i\a\ = 2{1 + X'^)), let us represent 
( |ii| ) in another form by reducing the integral to 
J^^^ arcsinh (A'^/^ cosx)dx and then using differentiation 
w.r.t. the parameter. We get for any A > 0: 



\Q/2 

4 / arctan x 



T^Q Jo 



lnA+ — 



X 

arctan x 



dx — 



(12) 



-dx. 



X 



In what follows, we always consider (the most inter- 
esting from the point of view of bandwidth) case A = 1 . 
It is not, however, difficult to generalize the argument 
below to the case of any A. For A = 1, the expression 
( p^ ) reads: 7(£ei) = iG/irQ. In these terms, the Thou- 
less conjecture says that the total bandwidth behaves as 
87(eeJ asymptotically for large Q. The following more 
general fact appears to hold and is supported by numer- 
ics. 

Let W{x) be the total length of the image of [0, x] un- 
der the inverse of the transform a = a{e) (recall that 
the whole spectrum is the image of [—4, 4], i.e., the total 
bandwidth is = VF(-4) -I- W(4)). Then for any co- 
prime P and Q (including P changing with Q such that 
P/Q is not small) as Q ^ cxj: 

2 Z"''^' 

W{a) - 47(0-) ^ — / K{x/4)dx, a e [-4,4]. 
ttQ Jo 

(13) 

The last equation is just another form of (^ij) for A = 1; 

K{k) = Jq^'^{1 — sin^ x)~^^'^dx is the complete elliptic 
integral of the first kind; the relation is our conjec- 
ture. 

It is easy to derive ( p^ ) in the case P = 1, where 
we can use the known from semiclassics |]l3| asymptotic 
expression for a{e). In this case, the main asymptotic 
contribution to Vl^(fT) comes from the bands in a small 
neighbourhood of e — 0. These bands have width of order 
l/QlnQ, while the others are exponentially narrow. For 
jej < 1, P = 1, and Q 00, 



3 



CT(e)-4cosh(£Q/4)cos( £^ln^- 



2ar8r(-+z-j- — 



(14) 



Hence, the asymptotic contribution of an individual band 
to W{a) is 



27r 



QlnQ 



arcsm 



|a|/4 



cosh(eg/4)' 



(15) 



and the number of such bands in an interval dt — d(sQ /4) 
is - InQ — . Therefore, the sum of all contributions 



W(a) - 4^ / arcsin ^^^dt = 
■kQ Jq cosht 

4^ /' 'i^(x/4)dx = 47(a), 

■^Q Jo 



(16) 



which completes the derivation of (^_3|) for P = 1. 

If we could show that for a general P, W{a) ^ cj{(t), 
where c is independent of cr, then the value of c could be 
found using a result of Last and Wilkinson |^ that for 
any coprime P and Q, I'^'l^OI^^ = Q^^ 1 where 

cr(ei) = 0, i — 1,...,Q. This can be written as 
|W^(0)| — 1/Q (right or left derivative) because, ob- 
viously, \W^ia)\ = J:?=lW'i^^)\~\ where a(e,) = a, 
i = 1, . . . ,Q. On the other hand, since K{0) — tt/2, we 
have |7^(0)| — 1/4(5, which implies that c = 4. 

Putting our results and conjectures together, we can 
roughly say that the logarithm of Onsager's partition 
function, the mean of the Lyapunov exponents, and 
the asymptotic bandwidth are basically the same object 
which is universal in that for a fixed value of a, it does 
not depend on P. 
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